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ABSTRACT 

We derive an analytic model for the redshift evolution of linear-bias, allowing for interactions and 
merging of the mass-tracers, by solving a second order differential equation based on linear perturbation 
theory and the Friedmann-Lemaitre solutions of the cosmological field equations. We then study the halo- 
mass dependence of the bias evolution, using the dark matter halo distribution in a ACDM simulation in 
order to calibrate the free parameters of the model. Finally, we compare our theoretical predictions with 
available observational data and find a good agreement. In particular, we find that the bias of optical 
QSO's evolve differently than those selected in X-rays and that their corresponding typical dark matter 
halo mass is ^ 10^'^ h^^ Mq and <; 5 x 10^'^ Mq, respectively. 

Subject headings: cosmology: theory - dark matter - galaxies: halos - galaxies: formation - large-scale 
structure of universe - methods: N-body simulations. 



1. INTRODUCTION 

The distribution of matter on large scales, based on dif- 
ferent extragalactic objects, can provide important con- 
straints on models of cosmic structure formation. How- 
ever, a serious problem that hampers such a straight for- 
ward approach is our limited knowledge of how luminous 
matter traces the underlying mass distribution. It is well 
known that the large scale clustering pattern of different 
extragalactic objects (galaxies, AGN, clusters, etc) trace 
the underlying dark matter distribution in a biased man- 
ner (Kaiser 1984; Bardeen et al. 1986). Such a biasing 
is assumed to be statistical in nature; with galaxies and 
clusters being identified as high peaks of an underlying, 
initially Gaussian, random density field. Furthermore, the 
biasing of galaxies with respect to the dark matter distri- 
bution was also found to be a necessary ingredient of Cold 
Dark Matter (CDM) models of galaxy formation in order 
to reproduce the observed galaxy distribution (eg. Davis 
et al. 1985; Cole & Kaiser 1989; Benson et al. 2000). 

Furthermore, the bias redshift evolution is very impor- 
tant in order to relate observations with models of struc- 
ture formation and it has been shown that the bias factor, 
6(z), is a monotonically increasing function of redshift. In- 
deed, Mo & White (1996) and Matarrese et al. (1997) have 
developed a model for the evolution of the correlation bias, 
defined as the ratio of the halo to the mass correlation 
function, the so called galaxy merging bias model. This 
model takes into account, via the Press-Schechter formal- 
ism, the collapse of different mass halos at the different 
epochs. In this framework and at high z's one finds, in 
the EdS universe, that b{z) — 1 oc (1 -I- z) (see also Sheth 
et al. 2001). According to Matarrese et al (1997), if the 
only halos that exist at any epoch are those that have 
just formed via a merging process of smaller halos then 
the bias evolution, in the high density peak limit, becomes 
b{z) oc (1 -I- z)^ (see also Bagla 1998). There are many 
studies of the merging bias model in the context of a A 



cosmology (e.g. Jing 1998; Sheth & Tormen 1999; Sheth, 
Mo & Tormen 2001), while the bias has been found to be 
a key ingredient of the halo model (e.g. Peacock & Smith 
2000; Seljak 2000; Ma & Fry 2000). 

The so-called galaxy conserving bias model (Nusser & 
Davis 1994; Fry 1996; Tegmark & Peebles 1998; Hui & Par- 
fey 2007) assumes only that the different mass-tracer (act- 
ing as "test particles" ) fluctuation field is proportional to 
that of the underlying mass and it predicts a bias evolution 
according to: b{z) = 1 + (bo — 1)(1 + z) for Jim = 1, where 
bo is the bias factor at the present time. Coles, Melott & 
Munshi (1999) have developed a bias model within the hi- 
erarchical clustering scenario which gives interesting scal- 
ing relations for the galaxy bias. Another approach was 
proposed by Basilakos & Plionis (2001; 2003), in which 
the linear bias evolution was described via the solution of 
a second order differential equation, derived using linear 
perturbation theory and the basic cosmological equations. 
Also, in this model the mass tracer population is assumed 
to be conserved in time. 

In this paper we extend the original Basilakos & Plio- 
nis linear and scale-independent bias evolution model to 
include the effects of an evolving mass-tracer population, 
possibly due to interactions and merging. We then derive 
the dependence of the linear-bias evolution on halo mass, 
with the help of N-body simulations, in the framework of 
the concordance ACDM cosmological model. Additionally 
we compare our model with observations. 

2. LINEAR BIAS EVOLUTION MODEL 

In this section we describe our general linear bias evo- 
lution model which is based on linear perturbation theory 
in the matter dominated epoch (eg. Peebles 1993), either 
assuming that the mass-tracer population is conserved in 
time or that it evolves according to a (1 -I- z)" law. 

2.1. Case 1: Conservation of Mass-tracer Population 
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In Basilakos & Plionis (2001; 2003) we have assumed 
that in the evolution of the linear bias the effects of 
non-linear gravity and hydrodynamics (merging, feedback 
mechanisms etc) can be ignored (eg. Fry 1996; Tegmark & 
Peebles 1998; Catelan et al. 1998). Then, using the linear 
perturbation theory we obtained a second order differen- 
tial equation which describes the evolution of the linear 
bias factor, 6, between the background matter and the 
mass-tracer fluctuation field: 

yd + 2{d + H6)y + A7:Gpn,6y^0, (1) 

where b ^ y + 1 and 6{t) oc D{t). Note that D{t) is the 
linear growing mode (scaled to unity at the present time) , 
useful expressions of which can be found for the A cosmol- 
ogy in Peebles et al. (1993): 

5nn,E{z) p (1 + x) 



D{z) 



(2) 



2 J. E^x) 
and for the quintessence models in Silveira & Waga (1994), 
Wang & Steinhardt (1998) and Basilakos (2003). The solu- 
tion of eq. IH) provided our bias evolution model (Basilakos 
& Plionis 2001). The solution for the different cosmologi- 
cal model enter through the different behavior of D{t) and 
H{t). 

In order to transform equation ([1]) from time to redshift, 
we utilized the following expressions: 

' (3) 



dz HoE{z){l + z) 
E{z) = p^{l + zf + 



1/2 



(4) 

with being the density parameter at the present time, 
which satisfies f^A = 1 — f^m- Taking into account the lat- 
ter transformations, the basic differential equation for the 
evolution of the linear bias parameter takes the following 
form: 

^-P(z)^+Q(z)y = (5) 
az'^ az 



with relevant factors, 
1 



P{z) 



1 dE{z) 



2 dD{z) 



1 + z E{z) dz D{z) dz 



and 



Q{z) = 



in^{i + z) 



(6) 



(7) 



2E'^{z) ■ 

The general solution as a function of redshift for flat 
cosmological models (Jlm-I-OA = 1) was found to be (Basi- 
lakos & Phonis 2001): 

h{z)^CiE{z)+C2E{z)I{z) + l , (8) 

with 



I{z) = 



(1 



E^{x) 



-dx 



(9) 



The integral of equation is elliptic and therefore its so- 
lution, in the redshift range [z,-|-oo), can be expressed as 
a hyper-geometric function: 

"13 7 njy 

6' 2' 6'~07 



I(z) oc (1 + z) 



-1/2, 



(10) 



,(1 + Z)3. 

Note that the first term in the right-hand side of eq.®, 
which is the dominant one, has an approximate redshift 
dependence ^ (1 -I- z)'^/^ while the second term has a de- 
pendence ~ (1 -f z). 



In Basilakos & Plionis (2001) we compared this bias 
evolution model with the halo merging model (eg. Mo & 
White 1996; Matarrese et al. 1997) for z < 3 and found 
a very good consistency, once we fitted the integration 
constants Ci , C2 by evaluating our model to two different 
epochs. Also in Basilakos & Plionis (2003) we compared 
our ACDM solution, evaluated at z = and z = 3 us- 
ing the Hubble Deep Field (HDF) results (Arnouts et al. 
2002), with the Matarrese et al (1997) model and again 
found a good consistency. 

However, the consistent comparison of our bias evolu- 
tion model with the available observational and numerical 
data is only after we normalize our model to two differ- 
ent epochs, given by the data. In other words there is no 
independent normalization of our bias evolution curves. 

2.2. Case 2: Evolving Mass-tracer Population 

The assumption used in the previous section, that the 
mass-tracer number density is conserved in time, should 
be a gross oversimplification, since in most galaxy forma- 
tion models it is expected that interactions and merging 
is very important at early epochs. 

We now drop this assumption by allowing a contribu- 
tion from the corresponding interactions among the mass 
tracers and thus an evolution of the mass-tracer popu- 
lation. We derive again the corresponding equation (1), 
starting from the continuity equation and introducing an 
additional time-dependent term, ^'(t), which we associate 
with the effects of interactions and merging of the mass 
tracers. We also make the same assumption, as in our 
original formulation, that the tracers and the underlying 
mass distribution share the same velocity field. Then: 

(5 -h Vu ~ and 5g + \Iu + ^{t) ~ , (11) 
from which we obtain 

5-5g = -^ . (12) 

Although we do not have a fundamental theory to model 
the time-dependent ^{t) function, according to the nota- 
tions of Simon (2005), it appears to depend on the tracer 
number density and its logarithmic derivative as well as 
on the tracer overdensity: ^'(i) oc ^'(n, (1 -I- 5g)dhin/ dt) 
(see eq. 10 of Simon 2005 and our appendix). 

Now, since we are dealing with linear biasing, we have 
5g — bS and using 6 = y -I- 1, we get that d{yS)/dt = — 
Differentiating twice we then have: y6 + 2yd + yS = — ^. 
Solving for y6, using the fact that yS = —y5 — VE* and 
utilizing the differential time-evolution equation of 5 (eg. 
Peebles 1993): 

d + 2H5 - ^TlGpraS = (13) 

we finally obtain: 

yS + 2{5 + H5)y + Ai^Gp^^Sy = -2H^ - * (14) 

which is the corresponding equation (1) for the case of 
interactions among the tracers. 

If we transform the latter equation from time to redshift 
the basic differential equation becomes: 



with 



g-P(.)g + Q(% = /(.) 

(1 + z)^'{z)E{z) - 2E{z)^{z) 



HoD{z) 



(15) 



(16) 
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where the prime denotes derivative with respect to red- 
shift. 

Therefore, if one is able to find a partial solution yp of 
eq. (fT5)) . then due to the fact that we know already the 
solutions [E{z) and E{z)I{z) in section 2.1] for the homo- 
geneous case, the general solution for the interactive case 
should satisfy the following formula: 

v{z)^CiE{z)+C2E{z)I{z)+yp{z) . (17) 

Using some basic elements of the differential equation the- 
ory, the corresponding partial solution becomes 

f{x)E\x)I{x) 



yp{z) = E{z) 



E{z)I{z) 



' f{x)E\x) 



Ax 



dx 



(18) 



/o (1 + ^)^ 

It is obvious that if the interaction among the tracers is 
negligible (^' ~ 0) then yp ~ 0, as it should. 

In order to proceed we need to define the functional 
form of the interaction term "^{z) which is not an easy task 
to do. Due to the absence of a physically well-motivated 
fundamental theory, we parametrize the interaction term 
utilizing a standard evolutionary form: 

^{z) = AH^{l + zY (19) 

implying that 

{l + zYE{z) 



J{z)=A{v-2)- 



D{z) 



(20) 



where A and v are positive parameters (to be determined 
from simulations see next section). This is a reasonable 
approach since it can be shown that indeed our "^{z) term 
has such a dependence (see appendix). From a mathe- 
matical point of view, we have that for v > 2 the bias 
evolution becomes stronger than in the case of no interac- 
tions, especially at high redshifts, which means that due to 
the merging processes the halos (of some particular mass) 
correspond to higher peaks of the underlying density field 
with respect to equal mass halos in the non-interacting 
case. On the other hand, the v < 2 case corresponds to 
the destruction of halos of a particular mass, which results 
into a lower-rate of bias evolution with respect to the non- 
interacting case. Now, for the limiting case with v = 2 we 
obtain yp = Q, implying no contribution of the interacting 
term to the bias evolution solution, as in the case with 
\l/(t) = 0, which can be interpreted as the case where the 
destruction and creation processes are counter-balanced. 

3. PARAMETRIZING THE BIAS EVOLUTION MODEL USING 
N-BODY SIMULATIONS 

Our analytical approach gives a family of bias curves 
with four unknown parameters (Ci,C2, A, v). In order to 
obtain the behavior of b{z) we need to somehow evaluate 
these constants for different halo-masses (since different 
halo-masses result in different 6(0) and different rate of 
bias evolution) and we do so by using a ACDM, dark mat- 
ter only, simulation. 

3.1. Simulation Data 

The simulation volume is of a 500 Mpc cube in 
which a random realization of the concordance ACDM 
(Om = 0.3, 17a = 0.7, h = 0.72 and erg = 0.9) power spec- 
trum was generated with 512'^ particles (Ragone-Figueroa 



& Plionis 2007). The simulation was performed with 
an updated version of the parallel Tree-SPH code GAD- 
GET2 (Springel 2005). The particle mass is rup > 7.7 x 
10^" h^^ Mq comparable to the mass of one single galaxy. 
The halos are defined using a FoF algorithm with a link- 
ing length I — 0.17(n)~^/'^, where (n) is the mean particle 
density. This linking length corresponds to an overdensity 
~ 330 at the present epoch (z = 0). 

We estimate the bias redshift evolution of the differ- 
ent dark mass halos, with respect to the underlying mat- 
ter distribution, by measuring their relative fluctuations 
in spheres of radius 8 Mpc, according to: 

KM,z)=^-^^, (21) 

where the subscripts h and m denote halos and the un- 
derlying mass, respectively. The values of as,hiM,z), for 
halos of mass M, are computed at different redshifts, z, 
by: 




(22) 



where N is the mean number of such halos in spheres of 8 
Mpc radius and the factor 1/7V is the expected Pois- 
sonian contribution to the value of cr| ,j . Similarly, we 
estimate at each redshift the value of the underling mass 
(Tg^m. In order to measure (t|j- we randomly place A^rand 
sphere centers in the simulation volume, such that the sum 
of their volumes is equal to 1/8 of the simulation volume 
(A'rand — 7500). This is to ensure that we are not oversam- 
pling the available volume, in which case we would have 
been multiply sampling the same halo or mass fluctua- 
tions. The uncertainties reflect the dispersion in the distri- 
bution of (t| j , taking 20 bootstrap resamplings of the halo 
sample. Note, that these errors do not take into account 
the uncertainties introduced by cosmic variance, although 
for the majority of the halo mass ranges investigated it 
should be negligible, due to the size of the simulation vol- 
ume. Nevertheless, we verified this by investigating the 
effects of cosmic variance on the value of trg for volumes 8 
times smaller than our simulation box and found them to 
be significantly smaller than the corresponding bootstrap 
errors. 

For our present analysis we use 10 different redshift 
snapshots, spanning the range: < z < 5 and five dif- 
ferent halo mass intervals: 

• Sample Si. 7.7 x 10^^ H-^Mq <M< IO^^H-'^Mq 

• Sample S2: 1 x IO^^H^^Mq < M < 3 x 10"/i-1Mq 

• Sample 5*3: 3 x IO^^/i-^Mq < M < 6 x IO^/i^^Mq 

• Sample 5*4: 6 x IO^^H-^Mq <M <1x IO^'^H'^Mq 

• Sample S5: 1 x 10^^h~^MQ < M < 4 x IO^'^H^^Mq 

In Figure 1 we present with dots the simulation b{M, z) 
results for the different subsamples (solid points Si, open 
points S2, crosses S3, solid squares 5*4 and open squares 
5*5). As it is expected, the biasing is a monotonically in- 
creasing function of redshift with its evolution being signif- 
icantly stronger for larger halo masses (eg. Kauffmann et 
al. 1999; Somerville et al 2001 and references therein). By 
solid curves we present our bias evolution models (case's 
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1 and 2 at the left and right panels, respectively), which 
have been fitted to the N-body results (see below). 




3 4 2 4 



z z 

Fig. 1. — The bias z-cvolution for different halo masses in a con- 
cordance ACDM simulation (points). Continuous lines represent 
our bias evolution model, fitted to the simulation results: {Case 1 
is shown in the left panel and Case 2 in the right panel). Note that 
the different continuous lines corresponds to different halo masses 
and thus to different values of the fitted model parameters, shown 
in Table 1. 



the zero-point bias, 6o Qj for the five subsamples, are pre- 
sented in Table 1, while we plot the results, as continuous 
lines, in the right panel of figure 1. 

Our Case 2 model fits extremely well the numerical sim- 
ulation results even at high redshifts, which implies that 
our modification of the test-particle bias solution, to take 
into account interactions among mass-tracers, is an ex- 
tremely good approximation to the more elaborate halo 
bias evolution solutions, based on the Press-Schechter for- 
malism. Note also that the derived values of the constants 
of integration, (Ci, C2), are very similar in both of our mod- 
els (Case 1 and Case 2). 




3.2. Fitting our Model b{z) 

Firstly we fit to the simulation data our (Case 1) model, 
where the mass tracers population is conserved in time. To 
do so we utilize a standard likelihood procedure and 
compare the measured bias evolution, b{M,z), with that 
expected from our theoretical approach (see eq. [5]). We 
define the likelihood estimatoi0 as: C{c) oc exp[— x^(c)/2] 
with: 



b'iM,z) - b'{c,z) 



(23) 



-1.5 
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Fig. 2. — The relation between the measured constants and the 
median of each halo mass interval (Ci and C2 in the right and left 



panels, respectively), 
given by eq. (24). 



The corresponding fits, in both panels, are 



where c is a vector containing the parameters that we want 
to fit and ai is the derived uncertainty (see previous sec- 
tion). Note, that in this case c — (Ci,C2). We sample the 
various parameters as follows: Ci £ [0, 25] and C2 G [—2, 4] 
in steps of 0.01. 

In Fig. 1 (left panel), we plot our derived bias evolution 
fits (lines) to the simulation data (points) to find that in- 
deed they compare very well for z^ 2.5 — 3. However, it 
appears that at larger redshifts the simulation halo bias 
evolves more rapidly than what our {Case 1) model pre- 
dicts, especially for the low mass halos (see Si, S2 and S3 
subsamples) , which of course should be attributed to the 
assumption of halo number density conservation. How- 
ever, it is evident that this model describes well the bias 
evolution for the high mass halos (see samples 5*4 and S5), 
since such high-mass halos form quite late in the evolution 
of structure formation processes. 

These results strongly indicate that the interactions 
among the relative low mass tracers at large redshifts play 
an important role in the evolution of halo biasing. We 
now perform our likelihood analysis for the Case 2 model, 
ie., using the corresponding vector: c = {Ci,C2, A,u). We 
sample the (A,!/) parameters as follows: A € [10""^, fO^^] 
in steps of 10"'^; and ly G [2,3] in steps of 0.02. The pro- 
cedure is applied only for the lower mass halos (see Si, S2 
and S3 subsamples). The resulting best fit parameters and 

^ Likelihoods are normalized to their maximum values. 

^ For z = we get feo = b(O) = 1 + Ci -I- C2l{0). Note that for the concordance model Q_ 



We would like to caution the reader that the fitted pa- 
rameters, tabulated in Table f, are based on a concor- 
dance ACDM simulation with a high as (= 0.9) normal- 
ization. We have investigated the sensitivity of these pa- 
rameters to the power-spectrum normalization by running 
a WMAP3 normalized (as = 0.75) simulation (but of a 
smaller box size: L — 250 h^^ Mpc) and found small dif- 
ferences and only for the Ci and 6(0) parameters. We find 
for the WMAP3 normalization case that the value of Ci 
and 6(0) are higher by ~ 4% and ~ 6%, respectively, with 
respect to those presented in Table 1. On the other hand, 
the choice of the value of Vl^ affects only the constant C2, 
which for the case of different flat cosmologies should be 
multiplied by a factor (ri,„/0.3)'^/^ (see Basilakos & Plionis 
2001). 

Now, in order to derive a generic halo-mass dependent 
formula of our bias evolution model, within the ACDM 
model, we correlate the constants Ci , C2 , with the median 
mass of the corresponding halo sub-samples and find the 
following fits (shown in figure 2): 



Ci,2(M) ~ ai,2 



M 



1013/i-iMp 



(24) 



where ai = 3.29 ± 0.21, /3i = 0.34 ± 0.07 and 02 = 
—0.36 ± 0.01, P2 — 0.32 ± 0.06 are the corresponding best 
values with their la uncertainties. 

Finally, we need to emphasize that with our current 
model, assuming also that the extragalactic mass tracers 
(clusters, galaxies, X-ray sources, etc.) are hosted by a 
halo of a given mass, it is straightforward to derive ana- 



1 - Qa = 0.3 we have /(O) ~ 9.567 
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lytically their bias evolution behavior. 




1 2 3 1 2 3 



z z 

Fig. 3. — Comparison of our b(z) model with different obser- 
vational data. Different line types represent different halo masses 
(see section 3.1). Left Panel: optical galaxies (open points), Lyman 
break galaxies (solid triangle), EROs (star), DRGs (open triangle) 
and 2dF radio sources (filled pentagon). Right Panel: optically se- 
lected quasars (open points and crosses), soft and hard X-ray point 
sources (open squares and solid diamonds). In the insert we plot, as 
solid points, the high-z SSRS DR5 QSOs and the same b{z) model 
that fits their lower redshift counterparts (ie.. A/ ~ 10^"^ A^©)- 

4. COMPARISON WITH OBSERVATIONS & OTHER BIAS 
MODELS 

We investigate the extent to which our bias-evolution 
model can fit recent observational data of different extra- 
galactic sources and which is the DM halo mass for which 
such a comparison is successful. In performing such a com- 
parison we implicitly assume that, on average, each class 
of extragalactic sources is related to DM haloes of some, 
specific, mass. Therefore, to the extent to which our model 
fits the data, we will be able to identify the (average) mass 
of the DM halo which the different extragalactic sources 
inhabit. We will also compare our model with other re- 
cent bias evolution models to show the extent to which 
they compare. 

In the left panel of figure 3 we compare our theoretical 
prediction with observations using optical galaxies (open 
circles) from the VIMOS VLT Deep Survey (Marinoni et 
al. 2005). It is quite evident that our model (solid line), 
corresponding to a DM halo mass of ~ 10^^ h~'^MQ, repre- 
sents extremely well the z-dependence of the optical galaxy 
bias. It is quite interesting that the same b(z) curve fits 
accurately the Lyman-break galaxies (solid triangle) at a 
redshift z > 3 (Steidel et al. 1998; Adelberger et al. 1998; 
Kashikawa et al. 2006). 

Roche et al. (2002) using the ELAIS N2 field and Geor- 
gakakis et al. (2005) based on the PHOENIX deep survey 
find that the correlation length of the extremely red galax- 
ies (EROs) lies in the interval ~ 10 — 17h~^ Mpc for a 
slope 7 = 1.8. Thus, the corresponding EROs bias value, 
derived using: 

/ \ 7/2 

b(z) = I — I D^+^z) with 7 = 1.8 & e = -1.2 , 

is 6 ~ 4.1 ± 0.9 at a redshift z ~ 1.2 (shown as a star in 
the left panel of Fig. 3). While, Foucaud et al. (2007) 
find for the distant red galaxies (DRGs) in the UKIDSS 
Ultra Deep Survey, a bias parameter of 6 = 4j;J g at a 
mean redshift of z ~ 1 (see open triangle). Our model 



bias-evolution curve that fits quite well both of these re- 
sults (EROs and DRGs) is the one corresponding to a 
DM halo mass of ~ 7.7 x 10^^h~^MQ (dot dashed Hue). 
For comparison we also show b(z) curves for two more 
halo masses: ~ 4.2 x lO^^h"^ Mq (long dashed line) and 
- 2 X IO^'^H-^Mq (sort long line). 

Finally, Maglioccheti et al. (2004) find 5 ~ 1.85 (see 
filled pentagon) for the 2dF radio galaxies at z < 0.15, 
corresponding to the model bias of a ~ 10^'*'h~^ Mq DM 
halo. 

In the right panel of figure 3 we compare the results 
of optical and X-ray selected AGNs with our b{z) mod- 
els. The bias evolution of optical quasars by Groom et al. 
(2005) and Myers et al. (2007) based on the 2dF QSOs 
(open circles) and the SDSS DR4 (crosses) survey, respec- 
tively, are well approximated by our 6(z) model for a DM 
halo of ^ lO^^h^^ Mq (short dashed line). It is interest- 
ing to mention that also Groom et al. (2005), utilizing 
the Sheth, Mo & Tormen (2001) theoretical prescription, 
found that the median DM halo mass, for the expected 
quasar bias evolution, is ~ 10^^h~^ Mq (Porciani, Maglioc- 
chetti, Norberg 2004; see also figure 7 of Negrello, Maglioc- 
chetti & de Zotti 2006; Hopkins et al. 2007). 
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Fig. 4. — The bias parameter as a function of the halo mass at the 
present epoch. The solid line corresponds to our predictions, while 
the dashed and the long dashed lines to those of the Jing (1998) and 
Seljak & Warren (2004) models, respectively. 

Recently, Shen et al. (2007) using the SDSS DR5 data 
found that the high redshift optical quasars (2.9 < z < 5.4) 
are more clustered than their z ~ 1.5 counterparts. The 
corresponding bias values (solid points in the insert plot) 
are fitted extremely well by the same model b{z) function 
as their lower redshift counterparts (ie., for a DM halo 
mass of - IQ^^h-^MQ. 

Finally, the relatively strong clustering results of X-ray 
selected AGNs, based on a variety of XMM and Ghandra 
surveys (eg. Gilh et al. 2005; Basilakos et al. 2004; 2005; 
Puccetti et al. 2006; Miyaji et al. 2007), correspond to 
bias results shown in the right panel of figure 3 as open 
squares (soft band) and solid diamonds (hard band). The 
model b{z) curves that fit these results correspond to halo 
masses M > 5 x 10^'^ Mq, strongly suggesting that X- 
ray and optically selected AGNs do not inhabit the same 
DM halos. 

From the previous comparisons it is evident that our 
bias evolution model (Case 2) describes well the bias be- 
havior of different extragalactic sources. 
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We now turn to compare our generalized test-particle 
bias model {Case 2) with some recent merging bias mod- 
els. As an example, in figure 4 we compare our solution to 
that of Jing (1998) and Seljak & Warren (2004) at ^ = 0. 
The corresponding functional b{M) forms appear to be 
quite similar, although the predictions of Jing (1998) and 
Seljak & Warren (2004), with respect to our model, are 
lower by ~ 20% and higher by ~ 15% in the low and high- 
halo mass range, respectively. It is worth pointing out that 
for DM halo masses close to those related to optical galax- 
ies, the Jing (1998) as well as the Seljak & Warren (2004) 
models predict a strongly anti-bias picture at the present 
epoch (&o ^ 0.68), which appears to be in disagreement 
with observational results which indicate that &o ~ 1 (eg. 
Verde et al. 2002; Lahav et al. 2002). 




1 2 3 4 5 



z 

Fig. 5. — Comparison of our bias model for a DM halo mass of 
lO-'^ Mq (solid line) with the corresponding Croom et al (2005) 
empirical form (dot-dashed line) and the Hopking et al. (2007) 
model (dashed line). The points represent the optical QSO bias 
from the 2dF (open points), the SDSS DR4 (crosses) and the SDSS 
DR5 (solid points) surveys, while at z ~ 0.2 we plot the Asiago- 
ESO/RASS QSO survey results (Grazian et al. 2004). 

Finally, in figure 5 we again plot the optical QSO bias, as 
a function of redshift (see also Fig. 3), together with some 



bias-evolution model predictions. Our generalized solution 
for a halo mass of M ~ 10^'^/i~^Mq [solid line] can be com- 
pared with the empirical formula [dot dashed line] derived 
by Croom et al. (2005; see their equation 15) and with 
the model of Hopkings et al. (2007; see their equation 8) 
[dashed line] . The comparison shows that the three curves 
fit well the recent QSO bias evolution {z < 2.5). At higher 
redshifts {z > 2.5) our model evolves significantly more 
than the other two models, and it appears to be the only 
model that fits the recent, observationally derived, high-.2 
QSO bias. 



5. CONCLUSIONS 

In this work we extend our original Basilakos & Plio- 
nis (2001) bias evolution model, based on linear perturba- 
tion theory and the Friedmann-Lemaitre solutions of the 
cosmological field equations, to take into account the red- 
shift evolution of an interaction term that could affect the 
mass-tracer number density. Parameterizing our analyt- 
ical model with N-body simulations of the concordance 
ACDM cosmological model (fim = 1 — f^A = 0.3), we 
investigate the halo-mass dependence of our linear bias 
evolution model, with or without the interaction term. 

Assuming that the cxtragalactic mass tracers are hosted 
by a dark matter halo of given mass, we can predict ana- 
lytically their present bias parameter and its z-evolution. 
Inversely, fitting our bias evolution model to the observed 
bias of different cxtragalactic mass-tracers, we can identify 
the typical mass of the DM halo within which they live. 

Comparing our theoretical predictions with the observed 
bias of different mass-tracers (galaxies, AGNs, DRGs, 
EROs, etc) we find a very good agreement with our bias 
evolution functional form and therefore we also identify the 
typical mass of the DM halo which they inhabit. This has 
allowed us to infer that X-ray and optically selected AGNs 
inhabit DM halos of different mass (with X-ray AGNs be- 
ing associated with higher DM halo masses). 
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APPENDIX 
REDSHIFT DEPENDENCE OF vj/ 

With the aid of the Simon (2005) nomenclature we derive the expected redshift dependence of our interaction term 
In Simon (2005) the related "sink" term, (j){z), is modeled as the variation of the tracer number density, nt^, given 
by their mass conservation equation: 

^ + iv(vn„,) = (Al) 
at a 

where a is the scale factor. Using the expression ritr — «.tr(l + Str), in which fitr is the mean tracer density (depending 
on redshift), then the relation between cj) and our ^' function is given by: 



1 



ntr{t) 



(bit) ~ {I + Str) 



dritr 



(A2) 



Now, we can attempt to explore the functional form of our ^! function, taking into account that the evolution of the 
galaxy luminosity function, $(i), can be parametrized (fc ~ 3) as: 

$(L,z)dL = (1 + z)'=$(L,0)dL . (A3) 

Therefore, the corresponding mean density as a function of redshift becomes: 

fe 



$(L)dL = nt,(0)(l + z) 
On the other hand, the time derivative of the mean density is written: 

dfitr dfitr dz 



= -Hokntr{Q){l + z)''-^{l + z)E{z) = -kHontr{z)E{z) 



dt 



dz dt 



Therefore we have that 



kHoil + Striz))E{z) 



(A4) 



(A5) 



(A6) 



_ntr{z) 

In order to proceed, we make the following assumption (similar to that of Simon 2005), namely a Taylor expansion in 
ntr{z) up to the second order: 

= A + Bntr{z) + Cn1^{z) + 
with A^B and C constants. Doing so, eg I A6I becomes: 

\^{z)\ ^ 



or 



|*(z)|^(l + z) 



nt,(0)(l + z)'= 
A 



+ B + Cnt,^ (0)(1 + z)'' + kHoE(z) + kHoE{z)Str(z) 



B 



.nt^(0)(l + z)2'= (1 + 2)^= 

where in the framework of the A cosmology, we have used that 

6tr oc {l + z)-^F 



CntriO) 



kHoE{z) kHaE{z)F{z) 



(l + z)'' (1 + z) 



fc+i 



(A7) 
(A8) 
(A9) 



iiii — 

3' ' 6 ' 17^(1 + z)3 



(AlO) 



(see Silveira & Waga 1994) with F{z) being the hypergeometric fmiction which is an increasing function of redshift. 
Therefore, for k ^ 3 and at the high redshift limit, we have: 

kHoE{z)' 



|*(z)|^(l + z)^ 



while for the Einstein de Sitter universe: 

\^{z)\ 



(1 + z)' 



CntriO) 

CfitriO)- 



(1 + Z)3 

kHo 



(All) 



(A12) 



(l + z)3/2_ 

It is evident that the second term of the right hand-side of eqs. lAllI and IA12I is small at large redshifts, which implies 
that indeed, as we assumed in our model of ^ (eq. 19 of main paper), \l/(z) oc (1 + z)^. We also show analytically that 
we should expect ~ 3, and indeed we find, fitting the simulation data, a similar value, ie., v ~ 2.5. 
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Table A1 

Fitted Case 2 bias evolution model parameters and their la uncertainties, for different halo-mass ranges. For 

THE S4 AND Ss samples, THE EFFECTS OF HALO INTERACTIONS IS NEGLIGIBLE AND THUS THE Case 1 MODEL IS VALID. 



Sample 


M/h- 




Ci 


C2 








V 


Si 


7.7 X 10^^ - 


1.0 X 10" 


2.12±0.09 


-0.24±0.01 


1.15 


0.91 ±0.01 2 ±0.08 X 10" 


-3 


2.56 ±0.11 


S2 


1.0 X 10^3 - 


3.0 X 10" 


3.98±0.10 


-0.43±0.01 


1.30 


0.98 ±0.02 5 ±0.20 x 10" 


-3 


2.62 ±0.14 


S3 


3.0 X IQis - 


6.0 X 10" 


5.48±0.85 


-0.57±0.09 


0.50 


1.16 ±0.04 6 ±0.30 X 10" 


-3 


2.54 ±0.15 


Sa 


6.0 X 10" - 


1.0 X lOi'^ 


6.78±2.80 


-0.68±0.30 


0.20 


1.37 ±0.10 






S5 


1.0 X 10" - 


4.0 X 10" 


9.09±3.00 


-0.89±0.50 


0.20 


1.67 ±0.10 







